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04 ; Abstract 

We study radiative decays of D*j(2317) and D s j (2460) using light-cone QCD sum rules. In par- 



ticular, we consider the decay modes D* sJ (2317) -» £>*7 and D s j(2460) -» D^7, D*j(2317)7 and 
evaluate the hadronic parameters in the transition amplitudes analyzing correlation functions of 

o 

scalar, pseudoscalar, vector and axial-vector quark currents. In the case of -D*j(2317) — > D*j 
we also consider determinations based on two different correlation functions in HQET. The decay 
widths turn out to be different than previous estimates obtained by other methods; the results 
favour the interpretation of D*j(2317) and D s j(2460) as ordinary cs mesons. 
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I. INTRODUCTION 



The observation of two narrow resonances with charm and strangeness, D*j(2317) in the 
Dvir invariant mass distribution 



mass distributions 



Q, fl fl fl B S Q and £> sJ (2460) in the D*n° and D sl 
[2, lil 4 3, I<j li| , has raised discussions about the nature of these states 



and their quark content 10|. The natural identification consists in considering these states 

as the scalar and axial vector cs mesons respectively denoted as D s0 and D' sl . In the heavy 

quark limit m c — > oo such states are expected to be degenerate in mass and to form a 

1 + 

doublet having sf = - , with se the angular momentum of the light degrees of freedom. In 
that interpretation the two mesons complete, together with D sl (2536) and D s2 (2573), the 
set of four states corresponding to the lowest lying P-wave cs states of the constituent quark 
model. A chiral symmetry between the negative and positive parity doublets D s — D* vs 
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would account for the equality of the hyperfme splitting 



D s0 — D' sl , suggested in ref. 
in the two doublets. 

However, estimates of the masses of these mesons based on potential quark models gen- 
erally produce larger values than the measured ones, implying that the two scalar and 
axial-vector cs D s0 and D' sl states should be heavy enough to decay to DK and D*K and 
should have a broad width. On this basis, other interpretations for D*j(2317) and D s j(2A60) 
have been proposed, for example that of molecular states |l3j . Unitarity effects in the scalar 
DK channel have also been considered Q|. 

Radiative transitions probe the structure of hadrons, and therefore they are suitable to 
understand the nature of D*j(2317) and L> s j(2460) distinguishing among different interpre- 
tations Q Their rates can be predicted by various methods and the predictions can be 
compared to the experimental measurements. In particular, it has been suggested that, in 
the molecular picture, the D s j(2460) — > D*j(2317)j decay should be driven by the D* — > 
transition and should occur at a different rate with respect to the rate for a quark- ant iquark 
meson decay 3|; such a suggestion has to be supported by explicit calculations in view of 
the experimental observations. 

The radiative decay widths £^(2317) -> D* sl and ZXj(2460) D { * ] ^, ^(2317)7 have 



been evaluated using the constituent quar 
(VMD) ansatz in the heavy quark limit 



and the Vector Meson Dominance 



£ model 

10, In this paper we use a different method, 



light-cone QCD sum rules, an approach exploited to analyze many aspects of the heavy and 
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|l7, 18], including radiative decays 19, 2^j| 



light quark system phenomenology (17|, [18|, including radiative decays |19J,|2J| (for a review 



and references see 2l|). We apply the method starting from the identification of D*j(2317) 
with D s0 and _D s j (2460) with D' sl and we discuss results and related uncertainties. In 
particular, in Section|n]we consider the decay mode D*j(2317) — > D* / y and describe in detail 
the calculation of the transition amplitude, the input quantities used in the analysis, the 
numerical results and the sources of uncertainties. In Section ITTT1 we carry out a calculation 
of the same transition amplitude in the infinite heavy quark limit, discussing the deviation 
from the case of finite mass which is sizeable in the case of charm. The radiative modes of 
-D^j (2460) are analyzed in Sections IIVIIVH where we find different results with respect to 
those obtained by other methods. In Section EH] we discuss the differences; in spite of them, 
considering the available experimental measurements, we conclude that the description of 
_D*j(2317) and D s j(2460) as qq states is favoured. 

II. ££,(2317) - £>; 7 

The amplitude of the El transition D s0 — > D*j: 

( 7 ( g , X)D* s (p, X')\D s0 (p + q)) = ed [(e* ■ rj*)(p ■ q) - (e* • p)(fj* ■ q)} , (2.1) 

with e(A) and ^(A') the photon and D* polarization vectors, respectively, and e the electric 
charge, involves the hadronic parameter d which has dimension mass -1 . According to the 
strategy of QCD sum rules, the calculation of this parameter starts from considering the 
QCD and the hadronic expressions of a suitable correlation function of quark currents. 
We consider the correlation function 

F lt (p,q)=i [ rf 4 xe^( 7 (g,A)|T[jt( x )j ( )]|0) (2.2) 



of the scalar J = cs and the vector J M = cj^s quark currents, and an external photon 
state of momentum q and helicity A. The correlation function can be expressed in terms of 
Lorentz invariant structures: 

F»(P, Q) = F [(p ■ e*)q, - (p ■ g)e*] + ... . (2.3) 

In order to compute F (or F^) in QCD, we carry out the light-cone expansion x 2 — > of 
the product of the two currents in (j2.2j) . This involves non-local matrix elements of quark 



operators between the vacuum and the photon state which can be expressed in terms of 
operator matrix elements of increasing twist. For example, contracting the charm-quark 
fields in eq. (|2.2|) we obtain 

= I 7^1 / <?x 6 _ 2 ^ ( 7 ( g> A)|a(z) 7 ^ + m c )a(0)|0) 

fc M < 7 (g,A)|a(a:)s(0)|0) 
- ik a (~/(q, A)|s(x)a Ma s(0)|0) + m c (~/(q, X)\s(x)^s(0)\0) 



d^k 
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(2.4) 



the expressions of the photon matrix elements in terms of distribution amplitudes are col- 
lected in the Appendix. This kind of contributions is depicted in fig^a). Moreover, the 
light-cone expansion involves higher-twist contributions related to three-particle quark-gluon 
matrix elements, as depicted in fig^b); the expressions of the relevant quark-gluon matrix 
elements can also be found in the Appendix. 




FIG. 1: Diagrams involving photon distribution amplitudes. The dashed lines represent the two 
currents in the correlation function (|2.2|) . In (a) two-particle contributions and in (b) three-particle 
quark-gluon contributions are shown. 



In addition to the contributions of the photon emission from the soft s quark, we must 
consider the perturbative photon coupling to the strange and charm quarks, fig. 2 (a) and 
(b). It produces an expression for F of the form: 

[ +C ° P P ( S ) 
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with 



= j - m , In ( S - < + m ; - *? m ]- m l ) + ^^Xi (s, ml,ml) 

3e s m c + m s As (s, m^, m^) /_ m 2 — m 2 



+ ^ 2 V ' i 1 — V^J + (s " c) (2 ' 6) 

(A the triangular function). Furthermore, nonperturbative effects when the photon is emitted 
from the heavy quark give rise to contributions proportional to the strange quark condensate, 
corresponding to the diagram in fig. 2 (c). 



(a) (b) ^ (c) X 

FIG. 2: Perturbative photon emission by the strange (a) and charm (b) quark. In (c) the strange 
quark condensate contribution is represented. 

The result is an expression of the correlation function (J2.2)) and of the function Fo in 
terms of quantities such as quark masses, condensates and photon distribution amplitudes. 
The sum rule for d is obtained by the equality of this QCD expression with a hadronic 
expression obtainted by a complete insertion of physical states. The two quark currents in 
(12. 2J) have non- vanishing matrix elements between the vacuum and D* and D s q: 

(0\4\ D *s) = fD*m Dt V, 

(D s0 \Jo\0) = f Da0 m Ds0 (2.7) 

so that can be written as 



(D 8 o|Jo|0)( 7 £>:|^o)(0|Jt|D; ) 

( m D s0 - (p + q) 2 )(™ 2 D* -p 2 



F^ = s „ ■; s 5 — 2 — ~ — oV + other resonances + continuum , (2.8) 



neglecting the widths of D* and D s0 . The sum rule follows after a double Borel trans- 
formation in — p 2 and — (p + q) 2 of both the QCD and the hadronic representation of the 
correlation function, that involves two Borel parameters, M 2 and Mf . The transformation 
allows to suppress the contribution of the continuum of states and of higher resonances, to 



5 



suppress the higher twist terms in the QCD expression of the correlation function and to 
remove all terms that are either independent of one of the two variables — p 2 or — (p + q) 2 
or depend on it only polynomially. The Borel parameters Mf and Mf are independent; 
we choose Mf = Mf since this allows, invoking global quark-hadron duality between the 
hadronic and the OPE expression of the correlation function above some threshold Sq, to 

2 2 

subtract the continuum in the QCD side through the substitution e~«^ — > e "iP — e ~ir^ in 
the leading twist term The masses of the charmed mesons involved in the transitions 
are close to each other, therefore the choice of equal Borel parameters is reasonable. The 
final expression of the sum rule for d is: 



D s0 Ds 
9 M% 



d 



m Ds0 f Da0 m Dt f D * 




8 p , v 

ds e p [sj 



(m c +m s ) 2 



+ e c e m 7 (ss) I 1 + 



m 



AM 2 2M 4 



+ e s (ss)(e _ M T — e~M I )M 



+ e s (ss)e 



- ^ (A(«o) - 8H 7 (u )) (l + 



M 2 



l—UQ 



dv I da g T{u§ — (1 — v)a g , 1 — u — vctg, ot gj 



+ I dv <ia g jF(ii — (1 — v)oi g , 1 — u — va g , a g 

n-uo Jo 



2e s f 3y m c e T&V v (u )) 



(2.9) 



where T = S - S - T t + T 4 - T 3 + T 2 + 2v{-S + T 3 - T 2 ), H^u) 



du'HJu'), HJu) 



du'h^v!) and ^ v {u) 



du'ip v (u') . All the distribution amplitudes are defined in the 



Appendix; u = - 2 - 

The sum rule (I2.9JI involves the meson masses, for which we use the experimental data, 
and the leptonic constants fo* and fz> a0 . For the former one, we put /d* = fjj a and use the 
central value of the experimental result fn B = 266 ± 32 MeV [22j. As for /d s0 , a sum rule 
obtained from the analysis of the correlation function 



U(p 2 )=t / d*xe ip ' x {0\T[J (0)4(x)}\0), 



(2.10) 



2 
n n 



f 2 

JD s0 



in 



D s0 



8tt 2 



so 



ds A 2 fs, m 2 m 2 . 



(m c +m s ) 2 



(m c + m s f 



e 



+ 



e M 



M 2 



M 4 



(sag s Gs) ml 
2 M 4 



'2.11 



allows to obtain fu s0 = 225 ± 25 MeV, using the parameters in the Appendix. 

From eq. ()2.9j) we can compute d varying the threshold So and considering the range of 
the external variable M 2 where the result is independent on it (stability region). In this 
region a hierarchy in the terms with increasing twist is observed, so that we can presume 
that the neglect of higher-twist contributions induces a small error. On the other hand, the 
perturbative term, which depends on both the light and the heavy quark charges, represents 
a sizeable contribution to the sum rule. 

In figElwe plot the curves corresponding to different values of sq- Considering the range 
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FIG. 3: The parameter d in the D s o — > D*7 decay amplitude ea . (|2.1j) versus the Borel parameter 
M 2 . The curves correspond to the thresholds sq = 2.45 2 GeV 2 (continuous line), so = 2.5 2 GeV 2 
(long-dashed line) and sq = 2.55 2 GeV 2 (dashed line). 



5 GeV 2 < M 2 < 7 GeV 2 , where the best stability in M 2 is found, together with the varia- 
tion of the threshold s , we get: —0.35 GeV -1 < d < —0.28 GeV^ 1 , corresponding to the 
radiative decay width 

T(D s0 -> D* s y) = (4 - 6) keV. (2.12) 

In (j2.12|) we have only considered the uncertainty in sq and M 2 , and we have used the 
central values of the QCD parameters collected in the Appendix. Actually, such parameters 
represent another source of uncertainty. In particular, an important input parameter is the 



7 



magnetic susceptibility of the quark condensate, Xi f° r which we use the value determined 
in ref.j3|: x = —(3.15 ± 0.3) GeV~ 2 . A different value x = — 4.4GeV~ 2 , previously used in 
other sum rule analyses, would produce a 40% larger value of \d\. 

The result (12.12)) shows that the radiative decay occurs at a typical rate for this kind of 
transitions (a few keV's). However, the rate is larger by a factor of 4-5 than that obtained 
using VMD and the infinite heavy quark limit, and by a factor of 2-3 larger than the 
estimates based on the constituent quark model. It is interesting to investigate the reason 
of the numerical differences; aiming at that, we estimate d by light-cone QCD sum rules in 
the heavy quark limit, using an approach based on the heavy quark effective theory. We 
discuss such a calculation in the next Section. 



III. L> sJ (2317) -> £>*7 IN THE HEAVY QUARK LIMIT 

In order to determine the hadronic parameter d in eq. ()2.1|) when m c — > oo, we consider 
two different correlation functions: 

if(w,?-«) = < f ^e^-^( 7 ( ? ,A)|T[jt(x)J (0)]|0) (3.1) 

and 

Ff\w,q-v) = % [ rf 4 xe^-^'(7(g,A)|T[jt(x)J d (0)]|0) . (3.2) 



The currents in (j3.1|) and (|3.2|) are effective currents constructed in terms of the strange 
quark field and of the effective field h v of the heavy quark (in our case the charm quark) 
with four-velocity v. The effective field h v is related to the heavy quark field Q in QCD 

1 + 2* n 

through the relation h v = e tmQV ' x — - — Q (for a review see [23]). The current J M = hyj^s 
has the quantum numbers of a vector meson. On the other hand, the currents Jo = h v s 
and Jd = h v (—i)'y IJ 'D-t tJL s have both the quantum numbers of a scalar meson, since D ttl = 
9tfiaD a = (flVta — VfiV a )D a , D being the covariant derivative. The latter current has been 
proposed as better suited for describing scalar heavy-light quark mesons in the heavy quark 



24j . therefore it is interesting to investigate how it behaves in sum rules for radiative 



limit 

decays. 

The sum rules for d are obtained from (|3.1|) and (|3.2|) using the same procedure followed 
in SecJnJ namely considering the light-cone expansion and the hadronic representation of the 
correlation functions, making a double Borel transform in the variables u> and u' = u — q ■ v 
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that involve two Borel parameters £1,2, choosing E\ = E?. = 2E and invoking quark-hadron 
duality above some threshold vq. From ()3.1|) we obtain 



4 A+A+ I 3m s e s 
— — e 2E 



FF^ 



1A) 



dv e E In 



v — {y 2 — m 2 s ) 2 



v + [y l — mj) 2 



(ss) 

+ e s ^-Ex0 7 («o) (1 - e" 



-in 



(ss) / A(u ) 



AE V 8 



HJhq] 



(3.3) 



On the other hand, from the correlation function (|3.2|) we get: 
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+ 3m s e, 
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dv e E In 
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1 — e B 



[m s + 1/) 



+ e s (ss) 



1 



E 2 (x0 7 (m o ) + %X0 7 (n o )) 



1-e B (l + ^r, 



+ e s (ss)- 



16 



(A(tt ) + itoA'(iio)) - H^(u ) 



(3.4) 



Notice that in ()3.3|) - ()3.4|) only photon emission from the light quark contributes. In the heavy 
quark limit the current-vacuum matrix elements are defined as follows: (0|J M |D*(f , X))h — 
Ffj^(X), (0\j \D s0 (v)) H = F + , (0\J d \D s0 (v)) H = F d + (the subscript H indicates that the 
states are normalized as used in HQET; FW and F^ have dimension mass 3//2 and mass 5//2 , 
respectively). Moreover, A and A + are mass parameters defined as A = m^. - m c , A + = 
f^D s0 — Trie (in the heavy quark limit). We use the numerical values: F 
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0.35 GeV 2, 
U4, 26, El. 



F + = 0.45 GeV5, F+ = 0.44 GeV 2 and A = 0.5 GeV, A+ = 0.86 GeV 
In figlU (left) we depict the result corresponding to eq. ()3.3|) . Considering the region where 
d^ is independent of the Borel parameter E: 1.2 GeV < E < 1.6 GeV, and the variation 
of the threshold u , we obtain -0.16GeV _1 < d^ < -0.13 GeV" 1 . Therefore, we obtain 
in the heavy quark limit a value compatible with the value obtained by VMD in the same 
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FIG. 4: The parameters obtained from eg . (j.S.Hj) (left) and d^ from eq, (|«-S.4|) (right) versus the 
Borel parameter E. The continuous, long-dashed and dashed lines correspond to the thresholds 
= 1.1, 1.2 and 1.3 GeV, respectively. 



limit: d ~ —0.15 GeV , finite quark mass effects are large, and enhance the D s $ — > D*^ 
amplitude by at least a factor of two. 

From the second sum rule eq. ()3.4|) . taking into account the dependence on the Borel 
parameter E for the continuum subtraction, we obtain a smaller result, see fig. 0] (right). 
This is due to a nearly complete cancellation between two different terms, the perturbative 
and the leading twist term, and therefore it critically depends on the input parameters of 
the QCD side of the sum rule, making the numerical result less reliable. 

IV. D sJ (2460) -> D sl 

Coming back to the case of finite charm quark mass, let us consider three radiative decay 
modes of D' sl , the transitions into a pseudoscalar D s , a vector D* and a scalar D s q meson 
with the emission of a photon. The calculation of the decay amplitudes is analogous to the 
one carried out in Section |nj therefore we present only the relevant formulae. 

The decay amplitude of D' sl — > D s j: 

( 7 (?, X)D a (p)\D' gl (p + q, A")) = egi [(s* ■ T})(p ■ q) - (e* ■ p)(rj ■ q)} (4.1) 

with r/(A") the D' sl polarization vector, involves the hadronic parameter gi that can be 
computed considering the correlation function 

T,(p,q) =ijd*x e^(j(q,\)\T[Jl(x)J^0m ■ (4.2) 
The quark currents are J5 = ci^s and = 07^755; can be expanded in Lorentz-invariant 
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structures: 



T,(p,q) = T [(e*.p) qfl -{p.q)s;]+.. 
The sum rule for g 1; obtained from the function T, reads: 



m 2 D , +m% s 



9i 



e 2aP (m c 4- m s ) 
m D'jD' a m 2 D J Ds 



+ e c e m^(ss) 
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dv / doi g Tx{uo — (1 — f)a g , 1 — uo — fO! g , a 5 
1— u >^o 



(4.3) 



+ 2e s fz 1 m c e ^V v (u ) 



(4.4) 



where JF X = S + S — T\ — T 2 + T 3 + T 4 + 2v{— S — T 3 + T 2 ) and the spectral function p p is: 

3e s 



2m s In 



s — m 2 + m 2 s — \z (s, m 2 , m 2 ) 
s — m 2 + m? s + A 2 (s, m 2 , m 2 ) 



+ (m c - m. 



ra -m -s ,i 



A2(s,m c ,m s ) 



c) 



(4.5) 



Eq. (|4.4j) involves parameters already used in previous Sections and the photon DA collected 
in the Appendix; it also involves the leptonic constant f& defined by the matrix element 



{0\tf\D' sl ) = f D > si m D , 1% , 



which can be obtained, starting from the two-point correlation function 



(4.6) 



IV(p 2 )=* J d 4 x ^ x {0\T[J^{0)jf{x 



(4.7) 
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from the sum rule: 



si 

e m 2 
m 2 D , 



ds A 2 (s, m 2 , m 2 s ) 



(m c +m a ) 2 



2 - 



m, + m 2 + 6m c m s 
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, , m 2 m s mlm 2 \ (scrg s Gs) ml 

{ss) m c - fn + 



2M 2 2M 4 
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(4i 



We get /d^ ~ / Ds0 . 

The calculation of gi produces the curves depicted in figJHl Considering the range 
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FIG. 5: The parameter g\ in the D' sl — > D*7 decay amplitude, ea. (|4.1j) . as a function of the Borel 
parameter M 2 . The curves refer to the threshold so = 2.5 2 GeV 2 (continuous), so = 2.55 2 GeV 2 
(long-dashed) and sq = 2.6 2 GeV 2 (dashed line). 



3 GeV 2 < M 2 < 6 GeV 2 , together with the variation of the threshold so, we obtain: 
-0.37 GeV" 1 <gi< -0.29 GeV" 1 , and therefore 

T(D' sl -> D sl ) = (19 - 29) keV. (4.9) 

As for D s0 — > .D*7, the result of light-cone sum rules for the width of D' sl — > D s ^y is larger 
than previous estimates. We shall discuss this point later on. 

V. D sJ (2460) -> D* sl 

The calculation of the dimensionless hadronic parameter g% appearing in the D' sl — > D*7 
transition amplitude: 

( 1 (q,X)D* s (p,X')\D > sl (p + q,X")) = te g 2 e a p (7T r 1 a f,^e^q T , (5.1) 
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with J7(A') and rj(\") the polarization vectors of D* and D' sl , is based on the analysis of the 
correlation function 



T, v {p,q) = i / Ae^( 7 feA)|T[j;(x)^(0)]|0) 



(5.2) 



expanded in Lorentz invariant structures 



T^(p, q) = T A E^ aT E*°q r + T B p„e v p aT ^e w q r 
+ T c (p + q) u e afM(TT p a e* a q T + ... . 
A sum rule for g 2 is obtained from T4: 

fj-2 = - — \ I ds e~M^ p p (s) + e c m c (ss) 



(5.3) 
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(5.4) 
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with JF 2 = S + S + T\ — T 2 — T 3 + T 4 and JF 3 = .4. + V. The perturbative spectral function 
p p reads: 

/{s) = |!i m , mcln C s -< + m ;-^( S ,^, m ;A + (s ^ c) (5 5) 

47r V s — + mj. + A2 (s, m^, m~) / 

The result is reported in fig El Considering the range 4 GeV 2 < M 2 < 6 GeV 2 and the 
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FIG. 6: The parameter 52 m the D' sl — > -D^7 amplitude ea. (|5.1() versus the Borel parameter M 2 . 
The continuous, long-dashed and dashed lines refer to sq = 2.5 2 GeV 2 , so = 2.55 2 GeV 2 and 
so = 2.6 2 GeV 2 , respectively. 

variation of the threshold sq, we get —0.18 < g<i < —0.13, i.e. 

T{D' sl ^ D* sJ ) = (0.6 - 1.1) keV. (5.6) 

The small value of #2 is due to large cancellations between the various contributions to the 
sum rule (|5.4|) : perturbative, twist two and higher twist contributions, as shown in figd 
In particular, the contribution proportional to turns out to be 50% of the contribution 
proportional to the magnetic susceptibility of the quark condensate. In the cancellation the 
detailed shapes of the distribution amplitudes and the numerical values of the parameters 
are of critical importance; this sensitivity induces us to consider the result for g 2 as less 
accurate than the results for the other channels. 

VI. D sJ (2460) -> D* i7 (2317) 7 

The last radiative decay mode we consider for D s j(2A60) is the Ml transition D' sl —>■ -D s o7 
which is governed by the amplitude 

("f(q,\)D s0 (p)\D' sl (p + q,\")} = ieg 3 e af3aT e* a V Pp°q r . (6.1) 
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0.8 r— 

0.6- 
0.4 i 
c 0.2- — 

01 - 
-0.2 
-0.4 i 



2 3 4 5 6 7 
M 2 (GeV 2 ) 

FIG. 7: Contributions to the sum rule (|5,4[) for g2- The continuous line corresponds to the per- 
turbative contribution in fig. 2 (a,b), the long-dashed line to the term proportional to the magnetic 
susceptibility of the quark condensate \i the long-short dashed line to the contribution propor- 
tional to /3 7 and the dashed line to the contribution corresponding to fig. 2 (c). The threshold is 
fixed to s = 2.55 2 GeV 2 . 



The parameter g% can be evaluated starting from the correlation function 



W lt (p,q)=i / *xe^( 7 (q,\)\T[4{x)j£(p)]\0) 



(6.2) 



written as 



(6.3) 



We work out the sum rule for #3: 



93 



m Da0 f Da0 m D 'f D ' 

J (m c +m,) 2 



so 



ds e'M^p (s) + e c e~M J (ss) I 1 + 



m s m c m 2 s m 2 c 



+ 



2M 2 8M 4 



+ e s (ss)(e m 2 — e «^)M 2 x0 7 (mo) 

+ e ^e a (sa)[--A(«o)(l + )] - -ye^^KK^ 



l-UQ 



+ e M 7 e s (ss) / / da g Ti{uQ — (1 — i> )a s , 1 — «o — UQ! g , a fi 



+ / / da. g T±{uo — (1 — u)a ff , 1 — u$ — va g , a g ) 

'l-uo Jo 



(6.4) 
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with F A = S + S + Ti + T 4 - T 2 - T 3 + 2u(-5 + T 3 - T 4 ) and 



P P (s) 



3e s j (m c + m s ) A| + ^ ^ fa - m 2 . + m\ - \\ (s, m% m 2 ) ' 



4tt 2 



- s 



c) 



s — m 2 c + ml + Aa (s, mj?, m 2 ) 



(6.5) 



Considering the range 4 GeV 2 < M 2 < 6 GeV 2 and varying the threshold Sq we get: 




FIG. 8: The parameter in the D' sl — ► .D s o7 amplitude, ea. (|6.1|) . versus the Borel parameter M 2 . 
The curves correspond to the same thresholds as in fig. [5] and El 



-0.35 GeV 1 < g 3 < —0.27 GeV 1 , corresponding to 

T(D' sl -> D s0 7) = (0.5 - 0.8) keV. 



(6.6) 



VII. DISCUSSION AND CONCLUSIONS 



As seen in the previous Sections, the radiative decay amplitudes of the charmed mesons 
considered here, when evaluated by light-cone QCD sum rules, are determined by two main 
contributions, the perturbative photon emission from the heavy and light quarks, and the 
contribution of the photon emission from the soft light quark. Other terms represent small 
corrections. In general, these two terms have different signs, and produce large cancellations; 
this allows to understand the role of QCD parameters like the magnetic susceptibility \- The 
delicate balancing of the two contributions determines the difference between the radiative 
widths of charged and neutral mesons. 



In Table |U we collect the LCSR results together with the results of other methods 
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161 ] . With the exception of D s j(2AQ0) — ► D*y, the rates of all the modes are larger 
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TABLE I: Radiative decay widths (in keV) of .0*^(2317) and L> sJ (2460) obtained by light-cone 
sum rules (LCSR). Vector Meson Dominance (VMD) and constituent quark model (QM) results 
are also reported. 



Initial state Final state LCSR VMD [10, ,16] QM [15] QM [11] 



£^(2317) 




4-6 


0.85 


1.9 


1.74 


D s j (2460) 


D sl 


19-29 


3.3 


6.2 


5.08 




D*sl 


0.6-1.1 


1.5 


5.5 


4.66 




££,(2317)7 


0.5-0.8 




0.012 


2.74 



than obtained by other approaches. In particular, T(D s j(2460) — > D s j) turns out to be 
considerably wider; notice that D s j(2460) — > D s j is the only radiative mode observed so 
far, as shown in Table ITT1 The peculiarity in £> s j(2460) — > D s ~f is that the perturbative 
contribution to the sum rule is the largest term, while in the other cases the leading twist 
term is the largest one in the theoretical side of the sum rules. 

TABLE II: Measurements and 90% CL upper limits of ratios of £^(2317) and £> sJ (2460) decay 
widths. 





Belle 


BaBar 


CLEO [2] 


r(r»:j(23i7)-z»; 7 ) 

r(D* J (2317)^D s 7rO) 


< 0.18 [3] 




< 0.059 


r(D sJ (2460)^D s7 ) 
r(D 3j (2460HDf7r") 


0.55 ±0.13 ±0.08 [3] 
0.38 ±0.11 ±0.04 [4] 


0.375 ± 0.054 ± 0.057 [9] 
0.274 ± 0.045 ± 0.020 [7] 


< 0.49 


r(D sJ (2460)^D* 7 ) 
r(D sJ (2460HD*7r u ) 


< 0.31 [3] 




< 0.16 


r(D sJ (2460)-*D; j (2317) 7 ) 




< 0.23 [9] 


< 0.58 


r(D sJ (2460)- > D>") 





In order to quantitatively understand the data in Table ILT1 one should precisely know the 
widths of the isospin violating transitions D s0 — > D s n° and D' sl — > -D*7r°. In the description 
of these transitions based on the mechanism of r] — rr mixing 

0, Q, Q, Q one should 

accurately determine the strong couplings D s oD s r] and D' s0 D*r] for finite heavy quark mass. 
Considering the results in Tables ITHTl these couplings should be larger than obtained in the 
heavy quark and SU(3) limit, an issue which deserves further investigation. 
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The obtained dominance of D' sl — > D s ^f with respect to other modes, in agreement 
with observation, induces us to consider our results as consistent with the interpretation 
of D;j(2317) and D sJ (2460) as ordinary cs mesons. The D sJ (2460) -> 1)^(2317)7 decay 
turns out to be suppressed; in the molecular interpretation it would be somehow enhanced. 
The observation of all the radiative decay modes with the predicted rates would of course 
reinforce our statement; in the meanwhile, we can reasonably conclude that invoking non- 
standard interpretations is not necessary. 

Acnowledgments. We thank T. Aliev for discussions and for collaboration in the early 
stage of this work. We acknowledge partial support from the EC Contract No. HPRN-CT- 
2002-00311 (EURIDICE). 
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APPENDIX A: PHOTON DISTRIBUTION AMPLITUDES 



For completeness, we collect in this Appendix the light-cone expansions of the photon 

matrix elements relevant for the calculation of the radiative decays of D s0 and D' sl . We also 

collect the expressions of the photon distribution amplitudes and the numerical values of the 

related parameters, as reported in t 28|. In all the expressions e(A) is the photon polarization 

s* ■ x 1 

vector and e M = e* — , g^ v = g^ iu -(q IJi x u + 4^); the variable u is defined as 



q ■ x 



q ■ x 



u — 1 — u; G^ v is the dual field G^ u = \s^ ua pG a ^ . We neglect quark mass corrections, that 
have not been worked out for all matrix elements. 



<7(?, A)|g(a)<vg(0)|0) = -iee q (qq) (s*q v - e^) / du e luq ' x X M U ) + 77^ A H 



x 
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ieeJ^-{x v Ea - XuE u ) [ du e luq ' x hJu) 
2qx J 



(7(9, A)|g(x)7 M g(0) 



(i{q, A) I g(z) 7^754(0) 



(7(9, \)\q(x)g s Gf, u (vx)q(0) 



(7(9, X)\q{x)g s G lxu (vx)i^ 5 q(0) 



{l{q, >)\q{x)g s G^(vx)-i a ^q(Q) 



(l{q, X)\q(x)g s G fMU (vx)i'j a q(0) 



0) = eej^ / du e^») 



1 r 1 

0) = — A ee q h 1 6^e™q a x p / 
4 Jo 

0) = -tee q (qq)(s;q u -sM J Va t e l ^ +va ^ x S(a t ) 

0) = -tee q (qq)(s;q u -sM J Va t e^ a ^ x S{a % ) 

0) = ee q f^q a (e;q u - elq,) J Va, t e^ +va ^ x A( ai ) 

0) = ee q f 3l q a (e;q„ - eM J Va t e^+^^V^) 
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(l(q, X)\q(x)a a/3 g s G tMl/ (vx)q(0)\0) = 



ee q {qq)- 



+ 



Sfigavqp - £fi9i3uq a - <-> v) 



e a g^q v - e a g vj3 q^ - (a <-> f$) 



I 



Va i e i( - a9+va ' l)q - x T2(ai) 



+ 



- Qpx a )(e*q u - e*g M ) 



(Al) 



oti = {a q , a q , a g } and / U(a,j) = / d« g / dctg / — a q — a q — a g ). The phot 

J Jo Jo Jo 

distribution amplitudes (DA) have the following expressions: 



on 



7 (m 
A(m 



/i 7 (m 



ip a (u 
V(a q , a q , a 



= 6uu (l + <p 2 cl (2m - 1)) 

= 40m 2 m 2 (3A; - k + + 1) + 8(C 2 + - 3C 2 ) [uu(2 + 13mm) 
+ 2m 3 (10 - 15m + 6m 2 ) hiM + 2m 3 (10 - 15m + 6m 2 ) \nu] 
= -10(1 + 2£; + )Cf (2m- 1) 

= 5 (3(2m - l) 2 - 1) + — (15^ - 5crf) (3 - 30(2m - l) 2 + 35(2m - l) 4 ) 



(1-(2m-1) 2 ) (5(2M-l) 2 -l)^(l + ^<-^ 



5A0u^(a q — a q )a q a q a 2 



3Q0a q a q a g 
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S(a q ,aq,a g ) = 30a 2 (k + k + )(l - a g ) + (Ci + Ci")(l _ a g)(l _ 2a g ) 



+ C2[3(a 5 - a q f - a g (l - a s )] 



5(a ? ,a 5 ,a 9 ) = -30a 2 (k - k + )(l - a g ) + (Ci - Ci + )(1 - a g )(l - 2a g 



+ C2[3(« 5 - a q ) 2 - a g (l - a g 



Ti(a q ,a q ,a g ) = -120(3( 2 + (£)(aq - a q )a q a q a g 

T 2 (a q , a q , a g ) = 30a 2 g (a q - a q ) [{k - k + ) + (Ci - Ci + )(1 - 2a 9 ) + (2(3 - 4a s )] 



%(a q ,a q ,a g ) = -120(3(2 - C 2 )( a q ~ a q )a q a q a g 



%(a q , a q , a g ) = 30a 2 g (a q - a q ) [{k + k + ) + (Ci + Ci + )(1 - 2a g ) + (2(3 - 4a £ 



(A2) 



2 , , y 

> w 7 



The parameters in the distribution amplitudes are: = —(0.0039 ± 0.0020) GeV 
3.8 ± 1.8, ljA = -2.1 ± 1.0 H; k = O^Ci = 0.4, (2 = 0.3, <p 2 = k + = d + = C 2 + = (at 
the renormalization scale /i = 1 GeV) [19j. The other parameters in the QCD sides of the 



sum rules, at the same renormalization scale, are: m c = 1.35 GeV, m s = 0.125 GeV 
(ss) = 0.8(qq) (q = u,d), (qq) = (—0.245 GeV) 3 and (qgaGq) = ml(qq) with wig = 0.8 
GeV 2 . Finally, for the magnetic susce ptib ility of the quark condensate x we use the value 
X = -(3.15 ± 0.3) GeV" 2 obtained in \&\. 
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